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Abstract: Quantum scattering theory in one dimension, stationary and 
time dependent, is reviewed. The basic operator formalism is provided and 
its relation to reflection and transmission amplitudes is spelled out. Symme- 
tries, the unitarity of the S matrix, and analytical properties are examined. 
Especial atention is paid to the scattering of wave packets, in particular to 
formulae for the transmittance, wave functions, dwell times, decay behaviour, 
short and long time behaviour. The simple real potential case is examined as 
well as complex potentials, potentials with different asymptotic limits at plus 
and minus infinity and time dependent potentials. Numerical, approximate 
methods, and transfer matrix techniques are also discussed. 


Resumen: Este trabajo es una revisión de la teoría cuántica de coli- 
siones en una dimensión, en sus versiones “estacionaria” y “dependiente del 
tiempo". Se proporciona el formalismo de operadores y su relación con las 
amplitudes de reflexión y transmisión. Se examina el efecto de las simetrías, 
la unitariedad de la *matriz S" y sus propiedes analiticas. Se presta especial 
atención a la colisión de paquetes de ondas, en particular se obtienen fórmulas 
pàra la transmitancia, funciones de onda, tiempos de permanencia, y para el 
comportamiento asintótico a tiempos cortos y grandes. Se examina el caso 
sencillo de potencial real así como potenciales complejos, potenciales con dife- 
rentes límites asintóticos en más y menos infinito y potenciales dependientes 
del tiempo. Se discuten también métodos aproximados y numéricos. 


lEste trabajo obtuvo el Premio de la Academia Canaria de Ciencias (Concurso del 
año 1996) correspondiente a trabajos de investigación o de revisión de temas de interés 
científico sobre materias corespondientes al campo de la Fisica. 
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1 Introducción 


La teoría cuántica de colisiones estudia choques entre partículas o entre 
partículas y campos externos. Las colisiones son procesos localizados en 
el espacio y en el tiempo. Es util distinguir cualitativamente tres etapas: 
antes, durante y después del intervalo temporal en el que la interacción es 
efectiva. Antes y después de la interacción el sistema se mueve "libremente" 
con respecto a la parte del potencial que se anula asintóticamente. La teoría 
también trata procesos de “decaimiento” en los que el sistema comienza en 
la zona de interacción y posteriormente se desplaza a la zona asintótica. 

En los textos de teoria de colisiones el caso sencillo de las colisiones 
elásticas se presenta de forma invariable en un espacio de coordenadas tri- 
dimensional (típicamente discutiendo campos de fuerzas centrales), y no se 
presta ninguna atención a las colisiones en una dimensión, Podemos encon- 
trar algunos ejemplos de colisiones en una dimensión en libros de mecánica 
cuántica. Pero estos tratamientos ignoran el lenguaje y las poderosas técnicas 
de la teoría de colisiones, tales como la formulación en términos de operado- 
res, y se concentran en soluciones directas de la ecuación de Schrodinger en 
casos resolubles. analíticamente (por ejemplo igualando la función de onda y 
su derivada en los bordes de la “barrera cuadrada” ) para obtener amplitudes 
de transmisión y reflexión. La teoría de colisiones en una dimensión es el 
marco adecuado para ir más allá de estos casos sencillos. 

Muchos sistemas físicos pueden describirse en una dimensión (1D): La 
aplicación de la aproximación de la masa efectiva en estructuras semicon- 
ductoras de capas lleva a una ecuación efectiva en una dimension [1]. Mul- 
titud de recientes publicaciones se dedican a estudiar aspectos particulares 


de las colisiones monodimensionales en estas “heteroestructuras”. Algunos 


fenómenos de superficie, y las reacciones químicas en ciertas condiciones pue- 
den también describirse mediante modelos en 1D [2]. La ventaja evidente de 
estos modelos es su simplicidad. Debido a ella son útiles desde un punto de 
vista pedagógico y también como herramientas de investigación. Permiten 
poner a prueba hipótesis, nuevas ideas, métodos de aproximación y teorías 
sin complicaciones innecesarias o costosas computacionalmente [3.4]. Por 
último, estos modelos se emplean con frecuencia para examinar cuestiones 
fundamentales de la mecánica cuántica. 

La teoría de colisiones en una dimensión espacial presenta interesantes 
peculiaridades que hasta el momento no han sido recogidas en ninguna mo- 
nografía, aunque existe un buen número de publicaciones dispersas sobre el 
tema. Este trabajo pretende cubrir parte de esta laguna. El asunto es sufi- 
cientemente amplio como para tener que efectuar una selección de materias. 
La que aquí se propone trata fundamentalmente de cubrir todos los aspectos 
que un usuario ocasional de la teoría puede requerir, e incluye una colección 
de fórmulas que no se encuentran en los textos de teoría de colisiones. Al- 
gunas no habían sido obtenidas con anterioridad y otras son sólo conocidas 
por un reducido número de expertos. En cuanto al nivel de rigor matemático 
hemos adoptado el “método del físico”. Esto significa que han de verificarse 
un número de hipótesis en las aplicaciones concretas, tales como la existencia 
de los operadores de Moller, la completitud asintótica de estados de colisión. 
o ciertas propiedades analíticas de elementos de matriz de la resolvente. Para 
los lectores con inclinaciones matemáticas senalaremos los puntos donde aún 
se requiere una clarificación rigurosa; esto ocurre especialmente en la sección 
dedicada a las colisiones dependientes del tiempo. Esta monografía se de- 
dica al problema directo que estudia la solución de la ecuación de Schrodinger 


con un potencial de interacción supuestamente conocido. El problema inverso 
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consiste en recuperar el potencial a partir de informacion parcial o total de la 
solucion. Por supuesto, los dos problemas estan relacionados, pero el inverso 
requiere unas técnicas y un lenguaje propios [5-11] que no discutiremos. 
Una importante sección se dedica a colisiones “dependientes del tiempo” 
y procesos de decaimiento. A pesar de la inherente dependencia temporal de 
un proceso de colisión los textos usuales prestan mas atención a la versión 
“estacionaria” de la teoría, que estudia las soluciones de la ecuación estacio- 
naria de Schrodinger. Esto se debe en parte a que el experimento de colisión 
convencional de haces moleculares (dirigido a obtener secciones eficaces) se 
realiza en condiciones cuasi estacionarias, y también porque los estados es- 
tacionarios forman una base para analizar el proceso real dependiente del 
tiempo. En muchos textos la evolución de paquetes de ondas se relega a 
justificar en el límite monoenergético la derivación de expresiones para la 
sección eficaz obtenidas mediante el formalismo estacionario, y a discutir 
brevemente tiempos de vida de resonancias. Sin embargo, en procesos de de- 
caimiento de estados inestables el parámetro relevante es el tiempo de vida 
en vez de la sección eficaz, y de hecho es cada vez más importante considerar 
explicítamente la dependencia temporal en casos más generales. Los expe- 
rimentos modernos de colisiones incorporan laseres con pulsos de brevisima 
duración (femtosegundos) que hacen posible observar el movimiento de los 
paquetes de ondas en tiempo real, y disenar estados iniciales particulares 
para lograr comportamientos dinámicos específicos. La teoría debe adaptar- 
se a estas nuevas necesidades prestando más atención a la caracterización 
temporal. De hecho, debemos estudiar el proceso completo de colisión, que 
incluye la etapa transitoria (durante la colisión), y no sólo los regímenes 
asintóticos. La bibliografía tradicional de teoría de colisiones destaca la co- 


nexión entre los dos regímenes asintóticos pero desprecia lo que ocurre en la 
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región intermedia de interacción por considerarla “no observable". 

Pueden plantearse muchas generalizaciones del caso sencillo de una partícula 
sin estructura que se mueve bajo el efecto de un potencial real y local en una 
dimensión. Nos limitaremos aquí a las colisiones, descritas en una dimensión 
matemática, de una partícula sin estructura con un blanco, la barrera de 
potencial. Hemos preferido una presentación más pedagógica que compacta: 
en lugar de describir una teoría general válida para potenciales arbitrarios 
pero notacionalmente farragosa, estudiaremos primero el caso más simple con 
detalle. Algunas de las complicaciones se discuten en secciones posteriores de- 
dicadas a potenciales con niveles asintóticos diferentes a derecha e izquierda, 
potenciales complejos no hermíticos e interacciones dependientes del tiempo. 
Una interpretación amplia de la expresión “colisiones en una dimensión" 
podría incluir colisiones inelásticas y reactivas de sistemas compuestos por 
partículas limitadas a un movimiento rectilineo pero no estudiaremos este 


caso. 


2 Premisas básicas y notación 


En un marco no relativista, el operador Hamiltoniano para una particula 
sencilla (sin estructura) en una dimensión puede escribirse como H = Ho-4- V, 
donde 


2 
H = P (1) 
2m 


es el operador de energía cinética en términos del operador momento Pop, y 
V es el potencial. Supondremos en general que V es un operador local, es 


decir en representación de coordenadas tiene como elementos de matriz 


(z|V|z') = é(z — 2)V (z), (2) 
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y que V(x) se anula con suficiente rapidez (el criterio se especificará más 
adelante) para valores grandes de |x|. [El subíndice “op” se usa en aque- 
llos operadores cuánticos que puedan confundirse con números ordinarios o 


funciones.] Los estados propios de pop y Ho son ondas planas con momento 
definido |p), 
(nip) — hein (3) 


Están normalizados a la delta de Dirac, 


(plp) = ó(p — p^). (4) 


Las relaciones de cierre (o resoluciones del operador unidad 1,,) pueden es- 
cribirse en representación de coordenadas o de momentos, 
lop = [ dels) f. dplp){pl: (5) 
3 Teoria formal: Operadores abstractos y 
parametrizados 


En esta sección y en la siguiente trasladamos a una dimensión la teoría formal 
que se describe en las monografías dedicadas a las colisiones en el espacio de 
tres dimensiones. Los libros de Taylor o Newton [12,13] pueden servir de 
apoyo para aclarar el origen de las fórmulas o aspectos de rigor matemático. 

El estado normalizado de la partícula en el espacio de Hilbert se denota 
como w(t). En el pasado y futuro remotos evoluciona esencialmente de forma 
libre, y tiende (en el sentido de un limite fuerte) a los estados asintóticos de 


entrada y salida, dent Y Psal respectivamente, 


v(t) ET Pent(t) D OO (6) 
P(t) 2 salt) , t>o. (7) 
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Los operadores de la teoría de colisiones se formulan en versiones abstractas o 
parametrizadas que no debemos confundir. (Esto ocurre con cierta frecuencia 
a causa de una notación deficiente). En esta sección se proporcionan las de- 
finiciones y fórmulas básicas. Los operadores de Moller conectan los estados 


asintóticos con el estado real y, y S conecta los dos estados asintóticos, 


W(t) = ON+den(t) (8) 
v(t) = N psalt) (9) 
Psalt) = Sdent(t) - (10) 


Los operadores abstractos se definen mediante limites (fuertes) que implican 


tiempos infinitos, 


Qs. = UON a ae (11) 
T = "yn. (12) 
S = Mn. (13) 


T y S son los operadores de transición y colisión respectivamente (la forma 
más habitual de nombrarlos es operador “te” y operador “ese” respectiva- 
mente). Estos límites pueden existir cuando los operadores actüan sobre 
estados del espacio de Hilbert de funciones de cuadrado integrable. Sin em- 
bargo es conveniente considerar una extensión que pueda aplicarse a ondas 
planas, que están fuera del espacio de Hilbert pero constituyen una base util 


para sus funciones. Definiremos primero los operadores parametrizados 


e 
— 
u 
— 
I 


1 + Go(z)Top(z) (14) 
V + VG(z)V, (15) 


m 
ne) 
~ 
R 
— 
ll 


donde z es una variable compleja con dimensiones de energía, y se han intro- 


ducido las resolventes de H y de Ho, G(z) = (z — H)! y Go(z) = (z - Ho)! 
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respectivamente. Las relaciones entre G(z), Go(z) Top(z) y V son 


G(z) 
Top(z)Go(z) 


Go(z) + Go(z)Top(z)Go(z) (16) 
VG(z) (17) 


Las resolventes son singulares en el eje real positivo (con un corte de rama) 
y en polos sobre el eje real negativo (estados ligados) del plano z. 

La relacion entre operadores abstractos y parametrizados se encuentra ac- 
tuando con (11) sobre una estado de cuadrado integrable. Puede introducirse 


un factor de convergencia sin cambiar el valor de la integral, 


0 > A 
Qu. = lime f dt ett eiH t/h ¿¡Hot/h (18) 
Q_ = Jim e | dic e n SPERM (19) 


Usando ahora una relación de cierre en momentos se obtiene 


Nlo) 
T*|9) 


J 4p (E, +10) Ip) (pl) (20) 
J 4p Top( Ep + i0)Ip) (plo). (21) 


A diferencia de los operadores abstractos, la acción de los operadores Q( E, + 
20) sobre las ondas planas está bien definida. En particular, pueden obtenerse 
de esta forma estados propios de H, 

1 
E, +10 — Ho 


Llamaremos al segundo término de la parte derecha de esta ecuación integral 


p^) = Q(E, + i0)|p) = |p) + Top(Ep + 20)|p). (22) 


“onda difundida” para distinguirlo de la “onda plana” del primer término. 
Para obtener una representación de coordenadas y evaluar el comportamiento 
asintótico se requieren los elementos de matriz de la resolvente, 


l A um +i|p||z—x'|/h 
= == a 23 
th; =o h|pl- ( ) 


Asintóticamente (|z| grande y z' fijo), 
I , IM silpzi/h y 
al —————l2) = EH e 24 
Las dos formas de acercamiento al eje real, desde abajo o desde arriba, en 
(23) implican diferentes condiciones de contorno. Para |p*) la onda difundida 
se forma con ondas salientes, y para |p~) con ondas entrantes. 


Podemos introducir las siguientes resoluciones de los operadores 2, T y 


S: 
0% = | dplp*)(p| (25) 
T = V J dp |p*)(p| (26) 
S = Jf dpdp'|p'\p|p*)(01. (27) 


Estrictamente los operadores en las ecuaciones (25-27) no son idénticos a los 
de las ecuaciones en (11-13) puesto que el dominio de los parametrizados 
incluye las ondas planas. Sin embargo, al actuar sobre estados del espacio de 
Hilbert son equivalentes, de manera que para evitar una notación complicada 
usaremos el mismo símbolo. 

Si el potencial V(x) es real la colisión conserva la norma, (¢ent|ent) = 


(v|v) = (dsailGsai). Esto significa que los operadores de Moller son isornétricos, 
NN, = 1. (28) 


Los estados de colisión v (t) con asíntotas entrantes y salientes se alejan del 
potencial en el pasado y futuro remotos por lo que son ortogonales a los es- 
tados ligados {|®,;)}, localizados en la zona de interacción, a tiempos (positi- 
vos o negativos) grandes. Como el solapamiento de dos estados cualesquiera 


(v|6;) = 0 es independente del tiempo, el espacio de estados ligados es de 
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hecho ortogonal al de estados de colisión, es decir, al rango de los operadores 
de Moller. Supondremos siempre que los rangos de los dos operadores de 
Moller son iguales (subespacio de estados de colisión), y que el espacio de 
Hilbert puede dividirse en el subespacio de estados ligados y el subespacio de 


estados de colisión. Esta propiedad se conoce como completitud asintótica, 
quem. Ne (29) 


y se satisface por la mayoría de potenciales que se encuentran en la práctica. 
En esta expresión el operador “deficiencia unitaria" A proyecta al espacio de 


estados ligados, 


A-X)0 (30) 


Como consecuencia de la isometría (28) la normalización de los estados |p*) 


coincide con la de las ondas planas, 
(p*|p^) = &p — p). (31) 


Pueden obtenerse varias resoluciones de la identidad usando conjuntos com- 


pletos de estados propios de H, 


lop = A+ | dplp*)(p*|. (32) 


Otro elemento importante de la teoria es la “matriz S”. En representacion 


de momentos tiene la forma 
Sp = ap- p) 7 2iró( Ep E Ey )(plToplEp d: 20)|p’) : (33) 


Usualmente se introduce una matriz S “en la capa de energia’ eliminando una 
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función delta de energía presente en (33) como un factor común. Usando? 


ap = p’) = Po, — Er) (34) 
donde óp» es la delta de Kronecker, 
] A pay - 
f des 
PE lo sl DAD, (35) 


y Spp = |plm~"6(E, — E’ )Spp', la matriz de colisión dos por dos en la capa de 
energía (p|S|p') se define mediante 


f 211m — > 
(plS|p') = bpp — “Tp PTE HiO) p=-tp’. (36) 


4 Simetrias 


4.1 Invarianza de inversión temporal 


El operador de inversión temporal 0 actua sobre los estados propios de la 


posición y del momento como 
0p) -|— p), lz) = |z}. (37) 
Es antilineal y antiunitario, 
0c|p) = c"0|v) (38) 
gm BE 1. (39) 


(c es cualquier constante compleja). Como @ es antilineal la acción del ope- 


rador adjunto 0! viene dada por 


(9. 0'w) = (06, vy . (40) 
“En general 
m 


&E,— Ej) = i 


[6(p— p') + é(p- p)]. 
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(Para evitar errores es conveniente usar la “notacion de los matemáticos" 
para los productos escalares en vez de la notación usual de Dirac.) 
Suponiendo que el potencial es real se cumple que [0, H] — 0, y por tanto 


Ó cambia el signo de los operadores de Moller, 


Como consecuencia. 
S — gisig, (42) 
y (en la capa de energía) 
(p|S|p') = (-rlS| — p) (43) 
(p|T^|p) = (—p'|T*|—p). (44) 


4,2. Paridad 


Muchos modelos de potencial son simétricos con respecto a una posición 
central zr = 0. El operador paridad II es unitario y cambia el signo de los 
vectores propios de posición y momento, II|z) = |—2), y II|p) = H|— p). Para 
Hamiltonianos invariantes frente a transformaciones de paridad (potenciales 


pares) |I, H] = 0 yslI; Q3 0. Borstanto, 


(p|T*|p') = (=p|T*| — p^) (45) 
(»|S Cali 2)- (46) 


p) 


5 Estados propios del Hamiltoniano H 


Los estados propios de H dados por las ecuaciones integrales de Lippmann- 
Schwinger (22) se comportan asintóticamente, de acuerdo con (24), como 


una combinación de ondas planas con momentos positivos y negativos. Los 
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factores que multiplican estas ondas planas son las amplitudes de rellexión v 


transmisión según la siguiente tabla (supondremos por ahora que p > 0) 


d J exp(ipr/h) + R'(p)exp(—ipr/h), si £~ —oc (47) 
hi/2 | T'(p) exp(ipz/h), si I^ OC, 

TIN T" (p) exp(—zpz/h), si T~ —oO (48 
hi/2 | exp(—zpr/h) + R'(p)exp(?pz/h), si z ~ oo. i 


Para p > 0 estas funciones son las representaciones de coordenadas de (z|pt) 
y (x|(—p)*), respectivamente, correspondientes a una onda plana entrante 
desde la izquierda, (r|p), o desde la derecha (z| — p). T(p) y R(p) son 
amplitudes de transmisión y reflexión para estos estados estacionarios. Si 
se formara un paquete de ondas de cuadrado integrable dominado por |p*) 
[14], este paquete sería localmente parecido a la onda plana |p) antes de la 
colisión. Después de la colisión se separaría en dos paquetes, uno reflejado y 
otro transmitido, con probabilidades |R|? y |T|?. 

Sin embargo para p « 0 los estados en (47) y (48) corresponden, respecti- 
vamente, a (z|p ^), con onda plana saliente (z|p), y (z| —p”), con onda plana 
saliente (z| — p). Como antes la expresión "saliente" se refiere a un paquete 
de ondas dominado por |p}. Este paquete de ondas sería localmente pare- 
cido a una onda plana sólo después de la colisión, de forma que la asintota 
entrante debe combinar ondas incidentes desde los dos lados de la barrera. 
Esta combinación puede ser dificil de conseguir en la práctica, lo cual no 
impide la utilidad de estos estados como funciones base, y en general en apli- 
caciones donde se controla o selecciona el estado final en lugar (o además) 
de preparse el estado inicial. 

Nótese que T(p « 0) no es una amplitud de transmisión estándar. Sin 
embargo, continúa analíticamente la amplitud de transmisión T(p > 0) al do- 


minio p < 0, y mantendremos el mismo nombre, “amplitud de transmisión”, 
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independientemente del signo de p. De acuerdo con nuestra convención los 
argumentos positivos de las amplitudes corresponden a estados |p*), mientras 
que argumentos negativos corresponden a estados |p"). 

El coeficiente de reflexión, R = |R(p)|?, es el cociente entre la magnitud 


del flujo reflejado, 7", y el flujo entrante j*"'; mientras que el coeficiente 


de transmisión, T = |T(p)|*, es el cociente entre los flujos transmitido 77 


3 


y entrante? De las anteriores definiciones se desprende que jf = j*"'(1 — 


|R(p)|?) y 
total. 


T(p)|? = 1—|R(p)|?, en concordancia con la conservación del flujo 


5.1 Relación entre los coeficientes de reflexión y trans- 
misión con los operadores básicos 


Para potenciales de soporte finito formados por secuencias de pozos y barreras 
“cuadrados” las amplitudes T y R se obtienen “empatando” la función de 
onda y su derivada en los extremos de cada pozo o barrera. Evidentemente, 
en casos más complicados es necesario un procedimiento más general que se 
describe en esta sección. 

Comparando el comportamiento asintótico (a |x| grande) de los estados 
en (47) y (48) con el comportamiento asintótico en (22), y usando (24), las 
amplitudes R(p) y T(p) pueden relacionarse con los elementos de matriz en 
la capa de energía del operador de transición. Detallaremos un caso a modo 
de ejemplo. Supongamos p > 0 y x — oo, 

(49) 


, 
LI 


(z|p)1, (50) 


/ dz'(x|Go( E, + i0)|z") ("| T. (E, + i0)]p) 


Oz mi etPrlh 2zmi 


i TOME USA "5 20)|p) ae 


~~ — 


3T también se llama en ocasiones “coeficiente de transmisión", y J “penetrabilidad” , 
“transmisividad”, “transparencia de la barrera" o “coeficiente de penetración”. 
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El resto de los casos puede tratarse de forma similar, en particular 


217m , 
T(p 21-— m (51) 
Hemos usado la invarianza de inversión temporal, que implica ye = m a 


para escribir T"(p) = T'(p) = T(p). Las otras amplitudes vienen dadas por 


2mir i 
R'(p) —2. z _— 


f 2m17 , si 
R' (p) E : (52) 


De (51) se obtienen las relaciones 


Pp T(p,  p real (53) 
R"(-p* = R™(p). (54) 


Usando (36), (51) y (52) la matriz S toma la forma 


_ (Tp) R(p) | 
S(p) = Rip) T(p) (55) 


p puede ser positivo o negativo pero de hecho uno de estos casos, por ejemplo 


p > 0 contiene toda la información de acuerdo con (53), 


S(—p) (56) 
S(—p). (57) 


N A 

x EL 

ms em 
as] "3 
uM 


Una diferencia importante entre las colisiones en la recta real, —oo < T < oo, 
y las colisiones radiales en un semieje, 0 < r « oc, es que en las primeras 
la matriz S es una matriz unitaria 2 x 2 mientras que en las segundas es un 
nümero complejo de módulo unidad. En un caso el espectro es degenerado y 


en el otro simple. 
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5.2 Soluciones de Jost y propiedades analíticas 


Ya se han descrito dos conjuntos de soluciones linealmente independientes 


[p*). Otro par de soluciones átiles viene dado por las soluciones de Jost, 
Jim ir) 1 
lim. Va(z)e'*7/^ — 1 (58) 
[Se entiende que lim, 0 Vi (z)'e7*?*/^ = ip/h y lim, , os Va(z)'e'?*/^ = —ip/h.] 
Estas condiciones de contorno iniciales explícitas facilitan su cálculo numérico 
y por tanto la obtención de las amplitudes, pero sobre todo permiten exa- 


minar las propiedades analíticas de las amplitudes de reflexión y transmisión 


en el plano complejo. Se relacionan con las soluciones en (47) y (48) por 


$i = T(p)'|p'e) (39) 
Ya = T(p)'|p =) (60) 


Las correspondientes ecuaciones integrales se obtienen por medio de las fun- 


ciones de Green para las condiciones de contorno (38), 
Vila) = e" = ja sin[p(z — z^)/^]V (z')v (z")dz' 
pala) = e — 7 f^ sin[p(z — a /MIV(s)vh(r)de. (61) 


Una prueba de la igualdad entre T" y T' se obtiene evaluando el Wronskiano 
W (1, 2) en ambos lados de la barrera, W (41, Y2) = —2ik/T! = —2ik/T" 
[8]. 

A partir del comportamiento asintótico de las soluciones de Jost se en- 
cuentran expresiones integrales para las amplitudes de transmisión y re- 
flexión, 

mi fo 
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¡CEPA (63) 


Eh Js 
mh" = = O e** V(z^)y(')dz' (64) 
T or m PS —ékz' Y ( LY y, / / Des 
TEN Eh fie V(z)va(z')dz (65) 


El estudio de las propiedades analiticas de las funciones y amplitudes en el 
plano complejo de momentos es un campo fascinante pero por falta de espacio 
nos limitamos a resumir las ideas y los resultados más importantes [15]. 
Las propiedades analíticas son muy ütiles para investigar el comportamiento 
en umbrales energéticos, las colisiones “resonantes”, efectos transitorios. o 
evolución a tiempos muy breves o muy grandes. Es fácil probar que las 
funciones (61) y los coeficientes (62-65) existen en el semiplano superior [8] 


para potenciales que satisfagan la condición 


/ dz(1 + 2?)|V(z)). (66) 


La posibilidad de continuarlas analiticamente al semiplano inferior depende 
de las propiedades específicas de la función potencial. El punto clave es el 
comportamiento de las exponenciales en (65). Estas exponeciales crecen en el 
semiplano inferior sin límite. Sin embargo, si el potencial es de soporte finito 
las integrales existen en todo el plano excepto en polos sobre el eje imaginario 
positivo (estados ligados) y polos en el semiplano inferior (resonancias). Si 
el potencial decae exponencialmente, compensa el crecimiento hasta cierto 
valor del momento imaginario, lo que permite establecer una banda de ana- 
liticidad en el semiplano inferior. El uso de coordenadas complejas permite 
también extender el dominio de analiticidad excepto posiblemente en el se- 
mieje imaginario negativo. 

Cuando T(p)^! = 0 en el semiplano superior el Wronskiano se hace cero. 


las dos soluciones de Jost son linealmente dependientes y disminuyen expo- 
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nencialmente en el infinito, es decir, corresponden a un estado ligado. Estos 
ceros deben ser simples [8]. 
Suponiendo que el potencial obedece la restricción (66) se tiene para 
Smp 2 0 [8] 
T(p) = 1 + O(1/Ipl). (67) 


Si T(p = 0) = 0, en la proximidad del origen, 


T(p) = bp + ofp), 6# 0 (68) 


para Smp > 0. 

La función compleja T(p)~! puede usarse exactamente como se usa la 
función de Jost de la onda parcial | = 0 en colisiones de potenciales esféricos 
para obtener el teorema de Levinson. Este teorema relaciona el nümero de 


estados ligados n con la fase de la amplitud de transmisión [16-18], 


zn —1 (no hay estado cuasi-ligado en p — 0) 


mn (estado cuasi-ligado en p = 0) 


D e) = l (69) 


El caso excepcional T(p = 0) Z 0 no corresponde a un estado normalizado, 


y se asocia con un estado cuasi-ligado. 


5.3  Unitariedad y sus consecuencias 


La unitariedad de la matriz de colisión S, SS! = StS = 1 proporciona dos 


relaciones: de los elementos diagonales se obtiene 
[T (p) + |R™“(p)/? = 1, (70) 
y de los no diagonales, 


T(p)[R'(p)]* + [T(p]' R'(p) 2 0, p real. (71) 
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Expresando las amplitudes de transmisión y reflexión en términos de módulo 
y fase, 

T(p = |T(p)e*r? (72) 

R"(p = |R'(p)e'*&), (73) 


la ecuación (71) permite relacionar las fases, 
207 + (2n +1)7 = Op + Opi, (74) 


donde n es un entero [19]. Resolviendo la ecuación matricial Sas = A,a se 


obtienen dos valores propios, 
Ac - THER” (75) 
A = THER). (76) 
Al ser S unitaria, estos valores propios tienen módulo unidad, como se com- 


prueba con (74), y pueden expresarse mediante “desfasajes propios", A, = 


ues = 0,1, 


bo 
ói 


[97 + arctan(|T|/| R[)]/2 (77) 
[®r + arctan(—|T]/| R])]/2 . (78) 


El producto de valores propios se relaciona con la fase de la amplitud de 
transmisión, det S = exp 2i®r. Los autovectores a, proporcionan dos estados 
linealmente independientes con la propiedad de que las amplitudes de ondas 
salientes son iguales a las amplitudes de ondas entrantes excepto por el factor 
de fase exp(276,). Los dos posibles cocientes con esta propiedad entre las 
amplitudes de ondas entrantes de derecha e izquierda son +(R'/R’)!/?. Si el 


otencial es par se cumple que R = R! = R’, y estas soluciones se convierten 
p pie q , 
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en pares e impares con ondas entrantes |p) + | — p). Esta circunstancia se ha 
aprovechado para elaborar un formalismo paralelo al que se emplea para las 


ondas parciales en tres dimensiones [20-24]. 


6 Dependencia temporal 


La dependencia temporal del estado (x|v(t)) puede expresarse mediante la 
cuadratura 

(clo) = f dp (zlp*)e 9^ (p+ |u(0)) (9) 
Toda la información sobre la evolución del paquete de ondas está en principio 
contenida en (79), pero es conveniente analizarla en términos de conceptos o 
parámetros sencillos. 

La colisión típica comienza con la preparación de un estado a un lado 
de la barrera con momento medio dirigido hacia la misma. Supondremos 
que el estado está a la izquierda y que tiene momento medio positivo. Tras 
el choque surgen dos paquetes de onda, uno transmitido y otro reflejado, 
con posiciones centrales retrasadas o adelantadas (segün las características 
del potencial y las energías implicadas) con respecto a un movimiento libre 
de referencia. La mecánica cuántica se manifiesta a través de una serie de 
fenómenos peculiares debidos a la deslocalización y coherencia de la función 


de onda: 


e En las colisiones resonantes parte del paquete de ondas permanece un 
tiempo mayor que lo habitual localizado en la zona de interacción. Si 
la resonancia corresponde a un pico en el coeficiente de transmisión y 
la anchura en momentos del paquete de ondas inicial es mayor que la 


anchura de la resonancia el paquete reflejado presenta dos máximos y 
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un minimo correlacionado con los momentos transmitidos preferencial- 
mente [25,26]. (Este tipo de correlación posición-momento se aprecia 
mejor empleando cuasi-distribuciones de probabilidad conjunta de po- 
siciones y momentos tales como la función de Wigner [27].) Si la reso- 
nancia corresponde, por el contrario, a un máximo en el coeficiente de 
reflexión, el mínimo aparece en el paquete transmitido. Las resonancias 
pueden en general asociarse a singularidades (polos) de las amplitudes 
de transmisión o reflexión en el semiplano inferior del plano complejo 


de momentos. 


El efecto tünel es otro de los efectos cuánticos paradigmáticos. Las 
potenciales aplicaciones del efecto tünel resonante, en el que la transmi- 
sión por debajo del umbral aumenta de forma significativa en torno a la 
energía de la resonancia, han impulsado todo un campo de investigación 
que estudia las propiedades y el diseno de estructuras semiconducto- 
ras. Un reto importante consiste en determinar las escalas de tiempo 
implicadas, pues la intuición clásica no es aplicable. Tan extrano como 
el efecto túnel, aunque poco comentado, es el efecto contrario según el 
cual paquetes con energía superior al máximo de la barrera son refleja- 
dos en parte. Esto ocurre con potenciales tan sencillos como un *pozo 


cuadrado”. 


Al discutir el efecto túnel en colisiones con paquetes de ondas hay que 
tener en cuenta que existe cierta ambiguedad en su definición. Como 
veremos en una sección posterior dedicada a la “transmitancia clásica”, 
la condición (E) < Vo (donde Vo es el máximo valor de V(z)) se con- 
sidera a veces como definitoria del efecto túnel. Pero esta condición 


no impide que parte de un colectivo clásico con la misma distribución 
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energética y valor medio que el estado cuántico atraviese la barrera 
con similar probabilidad que el cuántico. Para que el efecto tünel sea 
genuino es necesario que el paquete de ondas transmitido esté funda- 
mentalmente compuesto por momentos por debajo del umbral clásico 


(2mVo)#/?. 


e La colisión permite ademas el acceso a valores *prohibidos" clasicamente 
de determinadas variables, tales como el momento. En una colisión de 
un colectivo estadístico clásico con una barrera positiva que tienda a 
cero a derecha e izquierda se cumple necesariamente la desigualdad 


(para todo p, y todo t) 


f. do fo. < ie dp f(p, 0), (80) 


donde f(p,t) es la densidad de probabilidad de momento p a tiempo t. 


Sin embargo esta desigualdad puede violarse en el caso cuántico. 


6.1 ¿Qué es un estado asintótico? 


Puesto que los operadores de Moller son isométricos podemos escribir 


(p* (0) = (pIQ5 Q., |bene(0)) (81) 


en la ecuación (79). En muchas simulaciones numéricas el estado se coloca 
lejos del potencial en el instante en que se inicia el cálculo, convencionalmente 
i = 0. En un experimento, t = 0 puede identificarse con el tiempo de 
preparación. Al no existir solapamiento con el potencial en este instante 
la tentación es escribir (0) = ¢ent(0). En general esta substitución no 
es correcta. Recordemos el significado de $,,,. Es un estado’ que evoluciona 


libremente al cual tiende el estado real en el pasado remoto. Sin embargo, si el 
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estado inicial v(0) tiene componentes negativas parte del paquete se movera 
hacia el potencial al disminuir t, y finalmente chocará con él. Un verdadero 
estado asintótico localizado a la izquierda del potencial se caracteriza por 
tanto por tener sólo componentes de momento positivas. Paradójicamente 
esto es imposible desde un punto de vista estricto, ya que la localización en 
un semiespacio de coordenadas conlleva una deslocalización en momentos, 
pero en la práctica podemos acercarnos suficientemente a este caso ideal, por 
ejemplo mediante paquetes de ondas Gaussianos con varianzas seleccionadas 
de forma que el solapamiento inicial con momentos negativos y posiciones 


positivas sea despreciable. 


6.2 Algunas fórmulas útiles 


En una simulación numérica o en un experimento real el estado inicial a 
tiempo t = 0 no es necesariamente un estado asintótico, es decir no se puede 
identificar v(0) con Pen: (0) debido a la posible presencia de momentos nega- 
tivos. Supongamos únicamente que el estado inicial está bien localizado a la 


izquierda de la barrera de forma que se cumpla 


_ f (plv(0)) + (—plv(0))R'(p) p>0 
wq - (piv (0) T (=p) p «0 (82) 


Entonces la integral (79) puede separarse en dos partes, 
(ziv) = f dp(zlp*)e7 or^ (plu (0) (83) 


+ ji _ dpe IP (pl (0) [(zlp*)T(p) + (21 — 9*) (2) , 


donde se ha tenido en cuenta (53). La expresión entre paréntesis [...] es un 
estado propio de H. Usando (70) y (71), la consideración explícita de sus 


condiciones de contorno en +00 indica que este estado propio es |p~) (p < 0), 
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que continua analíticamente |p*) (p > 0) para momentos negativos. Como 


consecuencia, para valores grandes y positivos de z podemos escribir 


(alvet) ~ f dp (lpye EP (ppp (0))T (D). (84) 


y esto es todo lo que se requiere en muchas aplicaciones. La extensión del 
límite de la integral a —oo es un resultado notable y contrario a la intuición 
clásica. De hecho muchos autores truncan incorrectamente la integral en 


p= Uh 


6.3 La transmitancia 


La transmitancia T (4) es la probabilidad, para un estado dado y, de terminar 
con momento positivo en el futuro remoto. Excepto por el argumento usare- 
mos la misma notación para esta cantidad y para el coeficiente de transmisión 
(correspondiente a un momento fijo), 7 (p). 7 (V) y 7 (p) están íntimamente 
relacionados como se demuestra más abajo. [Los argumentos se omitirán 
cuando el contexto clarifique cuál de las dos cantidades se está manejando.] 
Como la distribución de momentos de un estado asintótico ¢,,; no cambia 


con el tiempo, la transmitancia puede escribirse como 


T) = f dp(pléca(0)) (és (Ip) (85) 


ION (I^) (86) 


Se obtiene una expresión sencilla suponiendo que a tiempo t = 0 el paquete de 
ondas no solapa con el potencial. En este caso se pueden introducir relaciones 
de cierre en el espacio de coordenadas y usar la expresión asintótica para el 


estado |p~) a la izquierda. El resultado es 


T(v) = [^ dp IOPI = [Carr ra), (87 
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donde se ha usado (53). 

La probabilidad de terminar con momento positivo es igual a la probabi- 
lidad de acabar en la parte derecha del potencial. El estado w se comporta 
finalmente como un estado libre sa, y puede demostrarse rigurosamente 


para el movimiento libre que cuando t — oo [28], 


Paté? = lim [^ dz leal (88) 
[dri = dim f^ dele (89) 


para cualquier coordenada finita a. 
Veamos otra interesante expresión para la transmitancia, como siempre 
suponiendo que el paquete está inicialmente a la izquierda de la barrera. El 


punto de partida es la ecuación de continuidad 


Ojv(z) _ 0J(z.t) 


E ub UE EL 
donde > 
dir z Im [V^ (z, t) Oy (z, t)/0x] (91) 
es la densidad de corriente cuántica. 
Integrando desde a hasta oo, y suponiendo J(oo,f) = 0, 
d NN, (a,t) m i 
v m J (at). (92) 
donde N,(a,t) es la norma a la derecha del punto a, 
N,(a,t) = ji dz |(z)f. (93) 


Finalmente, integrando de t = 0 a oo, y suponiendo N,(a,0) = 0 se obtiene 


T(V) = N,(a,00) = ji ^ dt J(a,t). (94) 
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6.4 Una medida de la duración de la colisión: El 
tiempo de permanencia 


La colisión de un paquete de ondas con una barrera de potencial queda des- 
crita completamente por la evolución del paquete de ondas desde el estado 
asintótico de entrada hasta el de salida. Sin embargo no siempre es necesa- 
ria toda la información contenida en vy(z,t). En muchos casos unas pocas 
cantidades características resumen toda la información relevante del choque. 
Una de estas cantidades es la transmitancia (que corresponde en colisiones 
de más de una dimensión a la sección eficaz). Otro elemento importante es 
la duración de la colisión. La medida estándar de la duración de la colisión 


es el tiempo de permanancia, o tiempo med:o de la colisión, 


ro(0,6:t,,t4) [^ dt [ dz li, (95) 


En principio las coordenadas a, b > a, y los instantes tı y tz > t, son arbitra- 
rios pero normalmente a y b se eligen para abarcar la región de interacción 
(de forma que V(z) sea cero o despreciable para x < a y b > a), tı se toma 
como tiempo cero (un tiempo inicial de preparación en el que el paquete de 
ondas aún no ha chocado contra la barrera), y t; = 00. 

El tiempo de permanencia se ha medido mediante la técnica de fotolumi- 
niscencia resuelta en el tiempo con un láser de picosegundos [29]; es además 
un parámetro importante para analizar el comportamiento de estructuras 
semiconductoras [30,31]. 

La interpretación de (95) como un tiempo medio de permanencia de la 
partícula en la región [a,b] (donde el promedio se refiere a los miembros del 
colectivo asociado con Y) no es obvia, puesto que no disponemos de trayec- 
torias en la interpretación ortodoxa de la teoría cuántica. Sin embargo, exis- 


ten varios argumentos formales que proporcionan (95) extendiendo al caso 
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cuántico el tiempo de permanencia clásico (por ejemplo mediante "integrales 
de caminos” de Fevnmann. trayectorias causales o de Bolin. o como el valor 
esperado de un operador tiempo de permanencia). No insistiremos aquí en 
estas cuestiones enralzadas en la deficientemente entendida dualidad onda- 
corpúsculo: independientemente de una hipotética interpretación estadística 
del tiempo de permanencia en términos de los miembros del colectivo. acepta- 
remos el tiempo de permanencia como una cantidad característica del propio 
colectivo (estado 4) que proporciona una definición razonable de la duración 
de la colisión, que al menos coincide con la expresión clásica. 

Veamos con más detalle la analogía entre (95) y el tiempo de permanencia 
clásico. Si un colectivo estadístico de partículas clásicas (que no interaccionan 
entre sí), descrito por una distribución de probabilidad en el espacio de fases 
f = f(p,q,t) normalizada a la unidad. colisiona con una barrera de potencial, 
el tiempo medio por partícula entre q = a y q = 6 puede escribirse como (las 
integrales de momento y posición van de —oo a oc, mientras que las integrales 


de tiempo van de 0 a oo si no se indica lo contrario): 


to = Jf Fao, po. 0) 7(a,b, qo. Po) dgodpo (96) 
= fdqo fapo [ f(ao po. 0) CPtla(ao. po. 1) —«] — Hla qo, po-1)— b)}ai 


2 ja n nig p.t) [H(q — a) — H(q — b)]dp = ju j o(q)dq 


donde 7(a, b, qo, po) es el tiempo entre a y b para la trayectoria con condiciones 
iniciales (qo, po), y H es la función de Heaviside, y (q. p) es la fase a tiempo ! 
de la trayectoria que comienza at (qo, po) a tiempo t = 0. De acuerdo con el 
teorema de Liouville dqdp — dqodpo y f(q.p,t) — f(qo. po, 0). En la ecuación 
(97), o(q) es la integral sobre momentos de la distribución f(q, p). es decir, 


la densidad de probabilidad. 
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Usando la notación Q = Q(a,b,t) = f? o(z)dz, el tiempo de permanencia 
se convierte en 
m f 9.5.04. (97) 
Integrando la ecuación de continuidad sobre z entre a y b, y sobre el tiempo 
entre 0 y t, Q toma la forma 


a(t) = [Ua t) - Jet) a= f F(a.bt di, (98) 


0 


donde J(x,t') es la densidad de corriente, F(a,b) = J(a) — J(b), y se ha 
supuesto que Q(0) = 0. Sustituyendo (98) en (97), 


TD 


fa [ 1 TUS fu fr — t) f (t) de 
lim “(= dt dim row- i is rie ae .(99) 


t'!—oo Jo t!— oo 


i 


Si además Q(t") decae más rápido que t"! [54] el tiempo de permanencia (97) 


toma la forma local (en [32] puede encontrarse una derivación alternativa) 
rp(a, b) = J [J(b, t’) — J(a,t’)] t' dt' . (100) 
Supongamos que a y 6 son tales que 
EIC dt! =T, "n J(a, t^) di! 1, DU J(a,1)d'=-R — (101) 


para al menos un instante t.. Físicamente el paquete de ondas pasa a través 
de a completamente antes del tiempo te, y depués de cierta duración con 
flujo despreciable la parte reflejada vuelve con una corriente negativa. En 


este caso (100) puede escribirse como 


ro = T (i) - (07 + RU (102) 


a 
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donde 


; IETLA MEE 


ir = —— 


J J(b,U) de 
(ye = f" Hatta (yz = 


(103] 


= Sd (ag td! 
J te, Ud! 


(101) 


En cada caso el instante medio de paso se obtiene normalizando la distri- 


bución de tiempos de paso. Teniendo en cuenta que R + 7 = 1, 


re = Qe O Te = (0% — (057 (105) 


a 


Siguiendo las referencias [33] y [34] se puede ver que las cantidades Tr y Tp 
tienden en el límite de distribuciones de momento muy estrechas a promedios 
de “tiempos de fase" definidos en la siguiente sección. Sin embargo, no 
coinciden en general con estos promedios, ni son necesarios estos paquetes de 


ondas especiales para su definición, sólo se requiere que exista un instante t.. 


6.5 Importancia de las fases. Retardos. Tiempos de 
llegada 


Si no hay estados ligados, puede obtenerse el potencial por inversión a partir 
de la matriz S o una de las amplitudes R' or R” en función del momento. 
Sin embargo, el conocimiento de los módulos, es decir de las probabilidades o 
coeficientes de transmisión y reflexión, no es suficiente. Las fases son también 
necesarias, y se asocian a propiededes observables con información sobre 
la dependencia temporal de la colisión. En esta sección usaremos ideas y 
resultados de las referencias [35-41]. 

Consideremos un paquete de ondas que incide desde la izquierda y choca 
contra una barrera (de 0 a d), y tomemos un intervalo espacial [a,b] más 


amplio que la barrera. A tiempo cero el estado es, en representación de 
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coordenadas, 


1 1/4 ; E 
HET Fea exp [ipox/h — (x — a.) /(46 J. (106) 
Está centrado en la posición rz, y momento p, = hk.. Su representación 


de momentos se denota como óo(p). La distribución de momentos inicial 
es f(p) = |¢o(p)|*, una distribución Gaussiana con varianza 0? = [h/(26))?. 


Lejos de la barrera, a la derecha, 


1 co 
Ur(z >> z5,t) = w. dp óo(p) T (p) e'"*—P9/^ . (107) 


y la izquierda, después de la colisión, 


pes E 
paT << tit v3! dp do(p) R'(p) e (?7*E9/^. (108) 


Hay un término adicional V;(r << zi,t), que describe la onda incidente. 


ent 


ent el efecto de Vg y su interferencia 


Antes de la colisión, para el cálculo de (t) 


con V; puede ignorarse. Después de la colisión, puede despreciarse V; para 


sal 
a? 


evaluar (1)?"', siempre que a sea una posición asintotica. 
Substituyendo las expresiones (106), (107) y (108) en los promedios tem- 


porales y haciendo uso de la delta de Dirac se obtiene 


O= x [Panes tr? rm quos 


xin d d 2 l gm Ls RT 
ty - x [" PIREZ f-a- a ean 


La cantidad 7F*(p) = m [b — z. + he, (p)] /p, se compone de un tiempo que 
emplearía una partícula clásica con masa m y momento f*en ir de ze to b, 


más un retardo de tiempo de fase mho'(p)/p. [De forma similar, el término 
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entre paréntesis en (110) es un tiempo que emplea la particula en viajar 
libremente de z. a z = a, con un cambio de signo instantáneo en r = 0, 
más un retardo. La integral da un promedio para el instante de llegada de 
la onda reflejada.] Por tanto, el tiempo de llegada (1);? puede verse como 
el promedio de los tiempos de fase r7^(p) sobre una pseudo probabilidad 
condicional de comenzar con momento p habiendo sido transmitido. (La 


evaluación del segundo momento proporciona términos adicionales sin inter- 


pretación semiclásica inmediata.) Usando 1; se obtiene 


2 m(a — Te) 


(yz = [^ aplés(p)i (111) 


La interpretación de estas cantidades como tiempos de llegada puede justifi- 
carse mediante modelos fenomenológicos de detectores (potenciales complejos 
absrbentes) [42]. 

La relación (109) es apropiada para examinar el *efecto Hartman" [43,44]. 
Hartman [43] estudió la evolución de un paquete de ondas con distribución de 
momentos centrada en p., chocando contra una barrera rectangular de altura 
Vo > p?/(2m), y anchura d. Encontró tres regiones que dependen del valor 
de d (véase también [45]). Para barreras anchas (opacas), los tiempos de fase 
asociados con k, por debajo de la barrera, tienden a una constante indepen- 
diente de d, 2m/(hkx), donde k = ,/2m(Vo — E). Cuando las componentes 
de momento por debajo de la barrera dominan la transmisión, el paquete de 
ondas transmitido parece atravesar la región del potencial en un intervalo 
de tiempo independiente de d. (“efecto Hartman"). Sin embargo, si el valor 
de d es muy grande (barrera ultra-opaca) las ondas planas con momentos 
por encima de la barrera dominan la transmisión, y resulta un comporta- 
miento clásico, es decir, el tiempo crece linealmente con d. Finalmente, para 


anchuras pequeñas, el tiempo de fase depende en general de d. 


4] 


Consideremos un paquete inicial Gaussiano (106) con el centro en z = z, 
y anchura espacial 6. Una estimación del valor de d que da la transición 
entre el efecto Hartman y el comportamiento cuasiclásico para obtenerse 
para cada valor de 6 igualando el factor |T(p)|* |óo(p)|? correspondiente a 
p = p. y p = pr, donde p, es el momento de la primera resonancia por 
encima de la barrera. Este procedimiento lleva a una curva que separa los 


comportamientos cuántico y cuasiclásico, 


5 PV -InIT()] | ñyred (112) 


|p- — pe] ae |p- T Pel’ 


donde ke = 4/2mW, — p2. Para 6 fijo, la transición es más abrupta a 6 
mayores como consecuencia de la distribución de momentos más estrecha. 
Un punto de vista complementario lo da la matriz de tiempo de retardo 


en la capa de energía propuesta por Smith [46], 


Q(B) = ise 02). (113) 


cuyos elementos diagonales son los tiempos de retardo para reflexión y trans- 
misión. Se llega a esta expresión substrayendo los tiempos de permanencia 
con y sin potencial [46,47]. La traza de (113) se relaciona con el cambio de la 
densidad de estados Ap(E) = Tr[6(E — H) — 6(E — Ho)], que es una cantidad 
fundamental para caracterizar el espectro continuo [48], de acuerdo con el 


"teorema espectral", 


Ap(E) = —- !S8mTr[G(E + i0) — Go(E + i0)] (114) 
E. m - dr(E) . 
= 7 Z aD =r E (115) 


Las versiones de ondas parciales y multicanal del teorema espectral se han 


probado rigurosamente [49] pero no conocemos ninguna derivación rigurosa 
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en una dimensión, aunque la justificación sigue presumiblemente líneas simi- 
lares. El resultado para una dimensión en términos de la fase de la amplitud 
de transmisión se ha formulado en [50] con argumentos plausibles que siguen 
a la referencia [51], y se ha confirmado con ejemplos explícitos [52]. La ültima 
expresión en la ecuación (115) se obtiene de la traza previa sobre Q usando 
las relaciones (70) y (74). Una derivación relativamente simple se obtiene 
siguiendo el razonamiento presentado en [53] para el caso tridimensional, y 


teniendo en cuenta los desfasajes apropriados para el caso monodimensional, 


(77) y (78). 


6.6 Comportamiento a tiempos grandes 


En esta sección consideraremos que el estado inicial v(z,0) puede o no so- 
lapar con el potencial a tiempo t = 0 como corresponde a un problema de 
decaimiento o a uno de colisión, respectivamente. Supondremos por simpli- 
cidad que no hay estados ligados, y t > 0 en todas las ecuaciones. 

Para discutir el comportamiento a tiempos grandes de la densidad de 


probabilidad se requiere la expresión del propagador en el plano complejo q, 


-i D i —iz 
(ze Pa) = 2 [ dle” (116) 
1 / 
w S3] is E. 


donde z = q?/2m y el contorno C va de —oo a +00 por encima de todas 
las singularidades de la resolvente. En ausencia de estados ligados C va por 
encima del eje real, donde se localiza el espectro continuo de H. Debido 
a la exponencial e7***/* en (116) el comportamiento cuando t — oo está 
dominado por la región en torno al origen, que es de hecho un punto de silla 


del camino de máxima pendiente de descenso de este factor exponencial. El 
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camino cruza el origen a lo largo de la diagonal de los cuadrantes segundo y 


cuarto. Introduciendo la variable auxiliar 


u=q/f, f—(1—i)y(mh/t), (118) 


la exponencial se transforma en e^", y u permanece real a lo largo del camino 
de máxima pendiente. 

El elemento de matriz de la resolvente (z|(z — H)^!|z^) que se define para 
Smq > 0 (o primera hoja de la energía) debe continuarse analíticamente en 
el semiplano inferior q (o segunda hoja del plano complejo z) para que este 
tipo de análisis pueda realizarse. Esto es válido en particular para potenciales 
de soporte finito. Si la función continuada analiticamente es analítica en el 


origen, posee una serie de Taylor 
(z|(z — H)^!|z^) = ao + a1q + azg? +... (119) 


con coeficientes a; que dependen de z y z'. Pero debido al factor (impar) q 
en (117), el primer término, ao, no contribuye a la integral (116). La fórmula 
asintótica para el propagador se origina por tanto en el segundo término y 


toma la forma 

3/2 

2 3 oe ae ee mh 

em x a= Ta E ) (120) 


Pueden darse contribuciones exponenciales debidas a los polos, pero no afec- 


(oen F^ z^ Es 


tan a este resultado asintótico. Cuando /(0) se hace cero aparece un com- 
portamiento asintótico t7? de la densidad de probabilidad. 
Mediante la descomposición de la resolvente de H 
1 1 1 1 
z— H z— Ho z — Ho z— H x 
1 1 m 


=i Ip? 
mc mue IR DE m 


(121) 


(122) 
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el operador de evolución U = e~*#*/* se separa en parte "libre" y "de difu- 


sion”, U = U, + Ua. 


(sls) = (alertan) = + f dget Na) (123) 
q l 
= hs 4 
4) = elii’) (124) 
(z|Ua|z’) = (z]e 8th — E $ dg «79^ 1g) (125) 
2:513 atm 
L4) = Sel uso) 7 12 (126) 


Pueden discutirse dos tipos de movimiento en una dimensión, movimiento en 
la recta real completa y en el semieje positivo. Mientras que en presencia de 
un potencial en los dos casos encontramos la dependencia t^? el movimiento 
libre es diferente. Consideremos primero el movimiento libre sobre la recta 
real. Usando la resolución de la identidad en términos de estados propios del 


momento e integrales de contorno se encuentra que 


n AAA, (127) 


1 
(alq le) = 
de manera que J;(0) = —1/h # 0, véase (124). Como consecuencia el compor- 
tamiento asintótico de la densidad de probabilidad para el movimiento libre 
en la recta real es t^!. Este es un caso importante en el que no se satisface la 


ecuación (119). Explícitamente, efectuando la integral en (116), se obtiene 


el conocido propagador 


gro emm. um 
f 1 


El análisis efectuado con varios potenciales modelo indica que la parte los 
términos libre y de difusión se cancelan, 1;,(0) = i/ = —1,(0) [54,55]. 
De diferente naturaleza es el movimiento libre restringido al semieje po- 


sitivo, r > 0. Un conjunto completo y ortogonal para este caso viene dado 
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por 


/ e 
(ips Iz sin(pr/h), r »0, p» 0, (129) 


que son las funciones propias de Ho, con valor propio p?/2m, normalizadas 
según fy (plr)(r|p')dr = 6(p — p'). La función de Green puede escribirse 


ahora como 


1 mt i(r+r’ ilr=r! 
as = E _ eilr—r'la/h) | (130) 


[Obsérvese la amplitud adicional en comparación con el resultado obtenido 


para movimiento libre sobre la recta real completa (127)]. I; es, según (130), 
L(a) = : [etre = esir=r"l0/8] (131) 


que se anula en el límite q — 0, lo que lleva a un comportamiento asintótico 
t^? de la densidad de probabilidad en contraste con la dependencia t^! del 
movimiento libre en la recta real. Explícitamente, el propagador es en este 


Caso 


LU Braet "X 
(r|U;|r’) = (=) let )?m/2th _ eitr c í (132) 


En los modelos analíticos examinados se cumple que /,(0) = I,(0) = 0 de 
forma que la densidad de probabilidad se comporta también asintóticamente 
como £^?. 

En resumen, en las colisiones en la recta real la combinación de contribu- 
ciones libre y de difusión cancela el término t7*, por lo cual la dependencia 
dominante de la densidad de probabilidad es t^?. En ningún caso el meca- 
nismo que explica esta dependencia es aisladamente el movimiento libre como 
se ha llegado a sugerir [56]. En el semieje r > 0 las dos componentes (libre 
y de difusión) tienen asintóticamente la dependencia dominante t^?. Podría 
darse la dependencia t^! si hubiera una singularidad en el origen. En el caso 


—oo < zr < oo un polo simple yace en el origen en ausencia de potencial 
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pero un potencial débil desplaza esta singularidad. Para valores arbitrarios 
del potencial sólo en casos excepcionales se encontrará el polo a energía cero. 
En colisiones de ondas parciales de tipo “s” las singularidades en el origen 


se conocen como resonancias de energía cero [12]. 


6.7 Comportamiento a tiempos cortos 


El decaimiento de estados inestables se ha estudiado de muchas maneras [57]. 
Una descripción ideal permitiría entender tanto el régimen exponencial do- 
minante como las desviaciones del mismo. Se ha progresado en esta dirección 
representando la amplitud de supervivencia A(t,w) = (v(0)|v(t)) como una 
suma discreta sobre términos resonantes [58,59]. 

Aqui discutimos el comportamiento a tiempos cortos del decaimiento de 
la probabilidad de supervivencia de estados cuánticos S — |A|?. Varios auto- 
res han descrito una dependencia a tiempos cortos t? de la probabilidad de 
decaimiento Pj. = 1 — S siempre que existan la energía media y el segundo 
momento, véanse por ejemplo los artículos relacionados con la *paradoja de 
Zenon cuántica" [60]. Se ha prestado menos atención al caso en el que estas 
condiciones no se satisfacen. Por otro lado Moshinsky y colaboradores sugie- 
ren una dependencia t!/? de la probabilidad de decaimiento a tiempos cortos 
[61-62]. Veamos cómo se resuelve la aparente contradicción. 

Cuendo el operador de evolución se expresa en términos de la resolvente, 


A(t, 5) toma la forma 


A AS —izt/h 5 i 
A(t, p) = (ple hy) = [dao gt) = z- [ dae M), 
(133) 


donde z = q?/(2m), y el contorno C va de —oo a +00 pasando por encima 
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de todas las singularidades de la resolvente. La función 


M(g) = tulg) (134) 


m 

puede evaluarse en el semiplano superior q, y luego ser continuada analitica- 
mente al semiplano inferior. Si la continuación existe, M(q) tiene en general 
un conjunto de singularidades de core, que dependen sólo del potencial, más 
posiblemente otras singularidades estructurales que dependen del estado par- 
ticular y. 

Es util deformar el contorno original de integración a la diagonal D de 
los cuadrantes segundo y cuarto del plano q para identificar las dependencias 
más importantes y por ser ventajoso computacionalmente. 


Supongamos que es posible un desarrollo de la forma 
ak 
M(g =>») ——~, Smg,«0 135 
(q) 2 ( = qk) qk ( ) 


(polos de orden superior pueden tratarse de forma similar). Aqui k = 
1,2,3--- son los indices de los polos. Al deformar el contorno C al con- 
torno D, los polos cruzados q, proporcionan contribuciones “resonantes” a 
A(t) que decaen exponencialmente con el tiempo (La colisión “resonante” se 


analiza con técnicas basadas en el plano complejo en [63]), 
Ex(t) = apet (mA) — qi eus (136) 


donde 
um qf, fm (1- i) (mh/t). (137) 


Independientemente de ser cruzados o no por la deformación de un contorno 
todos los polos contribuyen, debido a la integral a lo largo de una diagonal. 


Cada contribución se expresa en términos de la función wy+véase [64], como 


Dt) = -Z sign(Smu,) w|[sign(S*mu, Juz]. (138) 
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El término exponencial (136) puede añadirse e esta contribución para dar el 


resultado compacto [64], 
A(t) = » + Dx(t)] = 2 saxw(— —u;). (139) 


(Se sobreentiende que E,(t) = 0 para polos no cruzados al deformar el con- 
torno). En la primera forma de (139), el decaimiento exponencial en E; 
queda separado claramente de la “corrección” D,. Sin embargo la expresión 
compacta es muy util para estudiar el comportamiento a tiempos pequenos. 


La serie de Taylor de la función w [64] da una serie de potencias en t!/?, 


-ESÈ n (140) 
k ^ n=0 +1) 
Esto sugiere una dependencia t!/? de la probabilidad de decaimiento, como 
mantienen Moshinsky y colaboradores [61-62]. Por otro lado, la serie formal 
basada en desarrollar el operador de evolución, 
2 


— (VIB? lp) +---, — (141) 


A(t, 9) = (Ble!) = 1 — (IHW) — 302 
proporciona una dependencia t?, 
i? 
Pace = xa ((91H*|U) — (1 H10)%) +--+. (142) 


Sin embargo, los valores medios de H y/o potencias superiores de H pueden 
no existir. Son posibles varios comportamientos dependiendo de la exis- 
tencia de estos momentos. La cuestión de la posible existencia física de 
estados con momentos primero o segundo infinitos es un asunto controver- 
tido [65]. En principio, si son estados en el espacio de Hilbert, son estados 
válidos físicamente de acuerdo con la interpretación estándar del formalismo 
cuántico. Lamb propuso experimentos idealizados para crear estados arbitra- 


rios [66] mediante la creación de potenciales adecuados. Las actuales técnicas 
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de epitaxia de haces moleculares permiten el crecimiento de estructuras con 
formas de potencial muy flexibles [67], y proporcionan una via para acercarse 
a los experimentos ideales de Lamb. En cualquier caso, si los momentos no 
son estrictamente infinitos pero suficientemente grandes, el presente análisis 
seguiria siendo válido en cierta escala de tiempo. 

Consideremos las dos primeras derivadas de A a tiempo t = 0 primero a 
partir de (141), y suponiendo luego una dependencia de la forma A ~ 1+bf*", 


donde b y c son constantes finitas, 


dA =4 E 
"rhe T 6er], (143) 
dA? 1 B 

ames 7 cx) = be(e— 267)... (144) 


51 los momentos primero y/o segundo son infinitos los valores posibles de c 
están restringidos. En la siguiente tabla el comportamiento de la columna 
de la derecha es posible en principio cuando se satisface la condición de la 


columna izquierda, 


(V|H|v) 2 oo Pa. ^ t! (145) 
(p|H |b) < 00 and (V|H?|y) =00 Paec ~ t?” (146) 
(b|H |b) «oo and (V|H?|y) «oo  Paec t? (147) 


Pueden encontrarse ejemplos de comportamiento t!/?, 3/2 and t? a tiempos 
cortos de Pye, en [68]. Cuando la probabilidad de decaimiento se comporta 
como 1?/? or t? los coeficientes de potencias inferiores de t se anulan por las 


cancelaciones entre los diferentes polos. 
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7 Potenciales complejos 


Los potenciales complejos con parte imaginaria no nula aparecen de forma 
natural en colisiones *multicanal" al tratar de encontrar una ecuación efec- 
tiva para el canal de interés (típicamente el canal elástico) [12]. El canal 
seleccionado tendrá fuentes o sumideros físicos que pueden modelarse con un 
potencial complejo. Los potenciales complejos sirven también como modelos 
de detectores, y se usan en tratamientos fenomenológicos de sistemas disipa- 
tivos. Además, los potenciales complejos tienen una importante aplicación 
en propagaciones numéricas como absorbentes de paquetes de ondas en los 


bordes de la “caja” espacial finita en la que se confina al sistema. 


7.1 Formalismo 


En el caso hermitico, la evolución hacia adelante en el tiempo de (t1) a v(t;) 
viene dada por U (tz; tı) = exp [-1H (tz — t1)/A], (t2 > tı), y hacia atrás por 
U'(t2, tı) = exp [¿H(t, — t2). U es unitario y en una operación sucesiva de 
los dos operadores su efecto se cancela. Para un Hamiltoniano no hermitico 
U(t2;t,) no es unitario, ya que H no preserva la norma. Sin embargo se 
puede recobrar el estado inicial si el movimiento hacia atrás se realiza con 
el Hamiltonian dual H'. Definiendo la evolución hacia adelante bajo este 
Hamiltoniano mediante TIONI = exp [-iH'(t; — t,)/h], y hacia atrás por 
Üt (ta, t1) exp [iH (tz — t1)/h] encontramos 


O* (ta; ti)U(t2,t:) = 1. (148) 


Es importante tener presente para discusiones posteriores que si el movi- 
miento hacia adelante y hacia atrás corresponden a la acción física del Ha- 
miltoniano H, la exponencial para evolución hacia atrás en el tiempo se 


escribe en términos de Ht. Este punto se presta fácilmente a confusión. 
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Las anteriores consideraciones permiten generalizar al caso no hermitico : 
la unitariedad del operador de evolución (para Hamiltonianos hermíticos el 
símbolo” es irrelevante) e indican que podemos esperar también generaliza- 
ciones de las relaciónes usuales de la teoría de colisiones empleando H, Ht y 
sus funciones propias. 

La teoría de colisiones para Hamiltonianos no hermiticos podría formu- 
larse de acuerdo con diferentes convenciones de notación. Nosotros conserva- 
remos formalmente la estructura del operador S y la relación entre estados 
Ip?) y |p7) a través del operador de inversión temporal. En particular, el 
operador de Moller N, se define como en el caso hermítico, y admite la misma 


interpretación, pero (2_ viene dado por 


Q_ = lim e'F't/AgiHot/h (149) 


t—00 
asi que el movimiento hacia atrás en el tiempo que implica la segunda expo- 
nencial viene también determinado por H. De esta forma los estados |p~) se 


definen ahora mediante 

ip") - 0-Ip) - I) + Vo. (150) 
(Son estados propios de Ht.) El operador S que conecta el estado saliente 
con el entrante (de acuerdo con la evolución debida a H) viene dado de la 
forma habitual como S = Qin. La colisión no preserva la, norma y S no 
puede ser unitario. Sin embargo, un argumento similar al del comienzo de 
esta sección indica que la relación de unitariedad para 5 en el caso hermitico 
se generaliza introduciendo un operador de Moller que conecta los: estados 


asintóticos con el estado Y que evoluciona con Hf, 


Q, = lim in UTI: (151) 
à. = jim A gg eye (152) 
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Las relaciones (28) y (29) se convierten en 


HN, = lop (153) 
ont (154) 


I 
— 
o 
Y 
| 
-— 


Los estados |p*) tienen companeros biortogonales |p*) (se obtienen formal- 


mente a partir de las expresiones para |p*) intercambiando H por Ht y 


viceversa), 
1 
EIN d AD n 
Ip ^) Ip) + rap HR Ip) (155) 
| ] 
D) = _ _ _ 1 
IB 7) Ip) + I o e (156) 
a+ b 
(P lpt) = é(p'- p), (157) 


de forma que la unidad se expresa ahora en términos de una resolución bior- 


togonal, 
lop = JDA + f dplp*) (9*1 (158) 
El operador $ para movimiento bajo Ht es S = nus... y la relación que 


generaliza la unitariedad de S es 
1= Sts. (159) 


Las amplitudes de transmisión y reflexión para los estados propios de H, |p*) 
y |p) se escriben como es habitual (T, R), mientras que las de los estados 
propios de Ht se indican mediante T, R. 


Definiendo el operador parametrizado 


T(z) = V! 4 vt : E yt (160) 


a 


los resultados de la sección 5.1 pueden generalizarse fácilmente. En particular 


las ecuaciones para las amplitudes R y T siguen siendo válidas y se obtiene un 
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conjunto paralelo de relaciones mediante las substituciones R — R $9 
Usando el operador de inversión temporal se obtiene que 7" = T' = T, 
fe bares i y 


T(-p* = T(p) (161) 
R"(-py = Rp). (162) 


La relación de unitariedad generalizada (159) implica 


T(p)T (p)* + R"(p)R""(p)" 
R' (p) T(p) + T(p)' R'(p) 


1 (163) 
0, preal. (164) 


Una consecuencia inmediata es que el estado dependiente del tiempo 


(alo) = f dp (zlp*)e- 97^ (p+ (0), (165) 


viene dado por la misma expresión formal que en el caso hermitico. 

Un tratamiento diferente del presentado aquí, basado en la analogía con 
las colisiones de *ondas parciales" para potenciales pares puede encontrarse 
en [23]. 


7.2 Absorbentes 


Los potenciales absorbentes son modelos de detectores y sirven para elimi- 
nar los efectos espúreos de las paredes en cálculos numéricos. Por tanto es 
importante examinar algunas de sus propiedades. Supongamos que además 
del potencial complejo V4 existe otra interacción real V. Podemos definir un 
potencial absorbente como aquel para el que la norma de cualquier estado 


disminuye o se mantiene constante con el tiempo, 


E <0, paratodot. (166) 
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Es fácil comprobar mediante la ecuación de continuidad que una condición 
suficiente para satisfacer (166) es que la parte imaginaria del potencial Va 
sea negativa (o cero) para todo z. 

Un criterio alternativo se basa en la anulación de los coeficientes (y por 


tanto las amplitudes) de reflexión y transmisión para Va, 


R(p) 
T(p) 


0 (167) 
0 (168) 


I 


Una forma de asegurar (168) es imponer una barrera infinita en el punto 
zr = L. Si el potencial complejo tiene soporte finito, por ejemplo entre z = 0 
y z — L la condición (167) implica que la función de onda a la izquierda, 
xz <0, es igual a la función que existiría con el potencial V truncado en z = Q. 
Este resultado es significativo puesto que una de las propiedades peculiares 
de la evolución cuántica de paquetes de onda es que una superposición de 
ondas estacionarias con momento positivo puede dar lugar a flujos negativos 
a tiempos especiales [69]. Este efecto es en general pequeno o despreciable 
pero un Va estrictamente compatible con (167) debería dar cuenta de esta 
posible densidad de corriente negativa. Si se cumple (167) podemos igualar 
la derivada con respecto al tiempo de la norma a la izquierda con el flujo en 
$ 0, 
— M — Cee). (169) 
Esta expresión justifica operacionalmente (en términos de una medida con 
un detector ideal) la definición de tiempos de llegada discutidos en la sección 
6.5. 
Otras propiedades deseables para un potencial absorbente, sobre todo 
para las aplicaciones en cálculos numéricos, son (i) anchura del potencial 


L pequeña, y (ii) “robustez” frente a discretizaciones. La investigación en 
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este campo se dirige a encontrar potenciales que se acerquen tanto como 
sea posible a las anteriores condiciones mediante métodos variacionales o de 


inversion [70]. 


8 Niveles asintóticos diferentes 


Otra importante generalización del caso sencillo expuesto entre las seccio- 
nes 3 y 6, ocurre cuando V(z) tiende a valores asintóticos diferentes en los 
límites r — +00, V. y V, respectivamente. Esta asimetría es físicamente 
significativa en varias circunstancias. Por ejemplo, si se aplica una diferencia 
de potencial entre los extremos de la doble barrera de un diodo basado en 
el efecto tünel resonante, se crea un salto en el potencial. Este salto puede 
ajustarse para controlar la corriente que pasa a través de la estructura y lo- 
grar máximos o mínimos de corriente. También los modelos simplificados de 
reacciones químicas endo- o exo-térmicas presentan la asimetría correspon- 
diente a los distintos niveles de energía de los valles de reactivos y productos 
en las superficies de Born-Oppenheimer para el movimiento nuclear. 
Discutiremos la forma asintótica de las soluciones generalizando la sección 
5 y las relaciones para las amplitudes que se derivan de los Wronskianos. Para 


una energía E > max(V,, Vj) existen dos soluciones degeneradas linealmente 


independientes, 
.l fexp(ipz/h) + R'(p)exp(—ipiz/h), si z ~ —oo (170) 
hM? | T'(p) exp(tp,z/h), si Too, 
t T" (p,)exp(—ipız/ħ), si T~ —o0 (171) 
hi/2 | exp( —:p,z/h) + R’(p,)exp(tp,z/h), si z ^ ca. 


donde p, = [2m(E — V,)]/? y pı = [2m(E — V,)]?. Para valores positivos de 


Pr y pi, representan los estados |p?) y | — p?) respectivamente. Para valores 
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negativos, |p, ) y | - p; ). Estos pueden obtenerse de los anteriores mediante 
conjugación compleja de forma que las relaciones (53) siguen siendo válidas. 
Como el Wronskiano es constante encontramos 


i hh 


HEW) (el 97) = mü-IRP)-»ITP (172) 
wal pt) (eps) = -mIT'P--mü-IRP) (173) 
SS W((elpt),(2|— pf) = pil” = pT (174) 

PA W((elpt),(elpr)) = -nRT"-pR'T (175) 


Obsérvese que hay que distinguir entre amplitudes de derecha o izquierda. 
Sin embargo los coeficientes de transmisión de derecha e izquierda sí son 
iguales, 
flujo transmitido " 
T = ae = ae! (176) 
l 


flujo incidente 
flujo transmitido 
Toe = _ P. (177) 
flujo incidente Dr 
T! = T". Lo mismo ocurre para el coeficiente de reflexión, 


_ flujo reflejado 


RE = RP. (178) 


— flujo incidente 
Se define una matriz S unitaria mediante los cocientes entre coeficientes que 


multiplican a las ondas planas normalizadas a flujo unidad, 


TOJE Rp 
so) = | Rp) rent ) d 


Los dos ültimos resultados en (172) implican que las fases de las amplitudes 


de reflexión y transmisión se relacionan mediante 


Or: => Ör- = Or (180) 
Opi + Op = 11 +207. (181) 


DY 


9 Potenciales dependientes del tiempo 


Los potenciales dependientes del tiempo surgen como modelos sencillos en 
los que el efecto de uno o mas grados de libertad que no quieren conside- 
rarse explícitamente se simulan aproximadamente con una barrera V(z, t). 
También pueden corresponder a la acción de un campo oscilante acoplado al 
movimiento del sistema. En el caso general la evolución de un paquete de 
ondas en estas circunstancias puede calcularse con métodos descritos en la 
sección 10.1. La característica fundamental de estos choques es la aparición 
de energías que no estaban presentes en el paquete inicial debido a la ine- 
lasticidad. El caso en el que el potencial es oscilante admite un tratamiento 
especial [71]. Si el periodo es 27/w la ecuación de Schrödinger es invariante 


con respecto a la transformación 
t—t-4 20 : (182) 
Ww 
y aplicando el teorema de Floquet existen soluciones de la forma 
p(z, t) = e™*(z, t), (183) 
donde v es una constante, y ® tiene la periodicidad del potencial, 
$(z,t --27/w) = 9(z,t). (184) 
Desarrollando 9 en serie de Fourier, 


v(z,t) = = elt) p (4) (185) 


donde F, son los coeficientes de Fourier. Como 1h0/0t es el operador de 
energía, el término n-ésimo correspode a una energía hv + nhu, y (185) 


puede interpretarse como el estado correspondiente a una energía incidente 
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hv con componentes salientes, llamadas “bandas laterales", de energías dis- 
cretas debidas al potencial oscilante. Este formalismo se ha empleado para 
determinar el carácter sübito o adiabático de la colisión por efecto tünel. La 
frecuencia de oscilación constituye un parámetro de referencia con respecto 
al cual las colisiones se pueden clasificar como rápidas o lentas [72]. 

Una dependencia temporal de importantes consecuencias prácticas es el 
cambio de los niveles asintóticos que implica una variación de la diferencia 
de potencial aplicada a una estructura. En este caso es necesario determi- 
nar el tiempo de respuesta de una solución estacionaria de la configuración 
inicial hasta adaptarse a la configuración final. Este tiempo de respuesta 
limita la velocidad a la que idealmente podría trabajar un dispositivo semi- 
conductor. Existen publicaciones tratando este problema [73] aunque no se 
dispone de una teoría general. El problema en el que una barrera infinita, que 
actüa como compuerta, se elimina sübitamente se ha analizado con detalle en 
[74]. Azbel ha estudiado el efecto perturbativo de potenciales débiles depen- 
dientes suavemente del tiempo en las soluciones estacionarias del sistema no 
perturbado, en particular en la zona de efecto tünel, y describe resonancias 
(importantes amplificaciones de la transmisión) a energías y tiempos parti- 
culares [75]. Por ültimo, Martin y Sassoli de Bianchi han generalizado el 


teorema de Levinson para potenciales oscilantes [76]. 


10 Métodos numéricos y aproximados 


Los métodos aproximados son muy importantes en las colisiones que implican 
varias dimensiones espaciales debido a las dificultades para realizar cálculos 
exactos. En una dimensión el uso de métodos aproximados tales como W K B, 


la aproximación de Born [77] o métodos variacionales es posible en principio, 
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pero su interés está un tanto limitado por la relativa facilidad con la que 
se obtienen soluciones virtualmente exactas mediante métodos numéricos. 
No obstante las teorías aproximadas son ütiles cuando contribuyen a una 
explicación sencilla e intuitiva de los fenómenos o proporcionan resultados 


analiticos. 


10.1 Paquetes de onda dependientes del tiempo 


La mayoría de los métodos que se usan en la actualidad para propagar pa- 
quetes de ondas se han puesto a prueba primero en una dimensión [78], y 
existen varios programas disponibles. El método basado en el desarrollo de 
Chebychev del propagador combinado con la transformada rápida de Fourier 
es el más popular en dos o más dimensiones. Sin embargo, en una dimensión 
otros métodos son más rápidos y fáciles de programar [79], por ejemplo el 
método de Koonin [80]. Este método es simple, eficiente y estable. Se basa 
en el siguiente esquema unitario de propagación (A = 1), 
yoni, nm 
donde y” es la función de onda en el paso de tiempo n. El estado inicial 
más comun es el paquete de ondas Gaussiano (106). La posición espacial se 
discretiza mediante el índice 7, v(z, t) — v7. Los valores de z sé convierten 
en JAz, donde Az es la anchura de la celda espacial. De forma similar 
la variable tiempo toma los valores nAt, donde At es el paso elemental de 


tiempo. La derivada segunda espacial se aproxima como 
1 
i= Agi Pit — 2p; + Vj-1) + O(Az?). (187) 
El algoritmo se basa en escribir p"*? = xy — v^, donde 


(1+ iHAt/2)x = 29^ . (188) 


60 


La discretización espacial de la ultima ecuación lleva a un sistema lineal 
tridiagonal que puede resolverse mediante subrutinas adaptadas a este caso. 


El método se adapta con facilidad a potenciales complejos. 


10.2 Métodos de matrices de transferencia 


En un gran nümero de artículos se resuelve la ecuación de Schrodinger es- 
tacionaria mediante métodos de matrices de transferencia [81-90]. La idea 
básica es fragmentar el potencial en celdas. Las matrices de transferencia co- 
nectan la solución a ambos lados de cada celda. Este procedimiento se adapta 
muy bien a las series de pozos y/o barreras, que han despertado tanto in- 
terés por su conexión con la tecnología de semicondutores, dada la sencillez 
de las matrices implicadas. En realidad el método es aplicable en casos más 
complejos, ya que cualquier forma del potencial puede reducirse a pequerios 
fragmentos (de perfil plano o lineales) descritos por matrices sencillas. 

Kalotas y Lee [81] han desarollado una técnica basada en reemplazar el 
potencial por una secuencia de barreras o pozos de altura (interpolada) V; 
j =1,....,n y anchura w. Definiendo la matriz 


(189) 


—asinh aw cosh aw 


TNR coshaw  —(1/a)sinh aw ) 
jr] = 
donde a; = [2m(V; — E)/h*]!/? y usando para el producto de estas matrices 


la notación 


P = MaK [aw] (190) 


j=1 


la amplitud de transmisión viene dada por 


gnu 


E E. Po ME Bb 


(191) 
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Con técnicas similares Sprung, Wu y Martorell [82] expresan los coeficien- 
tes de transmisión y reflexión en términos de las amplitudes de las celdas 
individuales para potenciales compuestos por N celdas idénticas, lo cual es 
especialmente util en los potenciales conocidos como “super-redes”, que pre- 
sentan *mini-bandas" e interesantes propiedades de transporte como la po- 
sible resistividad diferencial negativa. 

Las matrices de transferencia también son ütiles para estudiar estacísti- 
camente colisiones en redes desordenadas (por ejemplo secuencias de barreras 
de anchura, posición o altura aleatoria) [83]. 

Rozman, Reineker y Tehver [84] derivan un sistema de ecuaciones di- 
ferenciales acopladas para las amplitudes de transmisión y reflexión de los 
potenciales “truncados” U*(zg) = U(zx)0(zg — x) de un potencial real arbi- 
trario U(z). Sean R'"(zo) y T(xo) las amplitudes de reflexión y transmisión 
del potencial truncado (en [84] se admite una posible asimetría en los niveles 
asintóticos). El sistema de ecuaciones viene dado en función de la variable 


a(x) = —imU(z)/ (kh?) por 


T'(zg) = a(xzo)T(zo)(1 + R'(zo)e?*7) (192) 
Rao) = a(zo)T'(zo)T" (zo)e?** (193) 
R''(zg) = alae "(IF R (ae y (194) 
con condiciones de contorno 
Ri(—00) = R'(-00)=0 (195) 
R(=) = 1. (196) 
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10.3 Colisiones múltiples 


El método de “colisiones multiples” [91-97] es similar al anterior. Consiste 


en descomponer el potencial en una suma de potenciales, 
yug (197) 


El operador de transición global Top puede escribirse como una serie en la 
que intervienen los operadores de transición individuales tu) intercalados por 
resolventes para el movimiento libre, 

Te” 5910 4- INDICAN +... (198) 

i ` i jži 

En casos favorables, esta técnica provee relaciones simples para las amplitu- 
des de transmisión y reflexión, y en general una imagen sencilla del proceso 
global en función de una secuencia de reflexiones y transmisiones en cada 
barrera. La amplitud total se obtiene como una suma coherente de las am- 
plitudes de los diferentes caminos posibles (clasificados según el número de 
interacciones), lo que lleva a una elegante visualización de efectos tales como 
las resonancias (que corresponden a energias para las cuales las fases de los 


múltiples procesos de colisión interfieren constructivamente). 


10.4 Ecuaciones diferenciales 


Es posible resolver directamente la ecuación estacionaria de Schrödinger me- 
diante subrutinas para sistemas de ecuaciones diferenciales de primer orden. 
El método se basa en descomponer la ecuación de Schrödinger de segundo 
orden en dos ecuaciones diferenciales acopladas de primer orden para las va- 
riables X = y(x) y Y 2 dy(z)/dz. Las condiciones de contorno se eligen de 


acuerdo con una de las funciones de Jost. 
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10.5 Aproximación clásica de la transmitancia 


Las transmitancias para paquetes de ondas dependientes del tiempo corres- 
pondientes a colectivos clásicos pueden ser muy próximas a las cuánticas. 
Este acuerdo puede darse también en la "zona de tuneleo" donde la energía 
promedio (E) es inferior al umbral clásico (o energía máxima de la barrera, 
Vo). Tal acuerdo se basa en la asociación estadística de la función de onda 
con un colectivo clasico, y se perdería si asociáramos una ünica partícula al 
paquete de ondas, como sugiere el teorema de Ehrenfest. El presente tra- 
tamiento supone un estado Gaussiano inicial con momento inicial medio p. 
[98]. 

Puede demostrarse que la transmitancia clásica depende sólo de dos pa- 


rámetros adimensionales f y r. 


a a 
r=-— I= mE) (199) 


donde m es la masa de la particula. También es importante, en el caso 
cuántico, el parámetro Up = (2mVo)/(h?a?), donde a^! es una anchura ca- 
racteristica de la barrera. La fórmula clásica para la transmitancia se obtiene 
asociando la Gaussiana inicial con un colectivo clásico con la misma distri- 
bución de momentos. La probabilidad de transmisión clásica es la fracción 


del colectivo con energía por encima del umbral clásico [98], 
T=1->P((1- M)/S] (200) 


donde P[(1— M)/S] es la integral de probabilidad de la distribución Gau- 
ssiana Z[(II — M)/S] con valor medio M = (fr)'? y varianza S? = (1— f)r. 
II es el momento adimensional II = p/(2mVo)!/?. 


De esta expresión puede concluirse que: 
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(1) La transmitancia no depende de la posición inicial del paquete. siem- 
pre que esté suficientemente alejado (que no solape inicialmente con el po- 
tencial). 

(2) La transmitancia no depende del parámetro Up ni de la forma del 
potencial. 


(3) La transmitancia depende sólo de los parámetros adimensionales f y 


(4) La transmitancia tiene un máximo en f — r, para f < 1. 

La validez del paralelo cuántico del primer punto se justifica observando 
las definiciones dadas en la sección dedicada a la transmitancia. Propaga- 
ciones numéricas de paquetes de ondas confirman la validez de las demás 
prediciones [98]. Se cumplen en el caso cuántico excepto en los siguientes 
supuestos: En el límite estacionario ( f próximo a 1), para valores pequenos 
de Uo, para potenciales discontinuos, y cuando las anchuras de las resonan- 
cia (si las hubiera) son mayores que la del paquete de ondas. Los resultados 
puramente clásicos son incluso mejores que métodos semiclásicos mucho más 
elaborados basados en la suposición de que el paquete evoluciona mante- 
niendo una forma Gaussiana [99]. 

Para obtener la fórmula de la transmitancia (200) no se ha supuesto más 
que el acuerdo con la distribución marginal inicial de momentos. El colectivo 
clásico puede especificarse imponiendo el acuerdo con la función de Wigner. 
La idea se ha explotado en una serie de aplicaciones en física molecular. Esta 
especificación es necesaria para obtener cantidades dinámicas y la evolución 


completa del colectivo clásico. 
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10.6 El método WKB y sus limitaciones 


E] método W K B proporciona una estimación sencilla del coeficiente de trans- 
misión. Para una barrera con dos puntos de retorno a « b a energía E el 


coeficiente de transmisión W K B viene dado por 
2 UE) 1/2 
Twp = exis |, del2m(V — E) (201) 


La mayor parte de los textos de mecánica cuántica advierten del fallo de esta 
expresión para energías muy cercanas al máximo del potencial y para po- 
tenciales que cambian con demasiada rapidez, y nos llevan a creer que para 
barreras opacas (anchas y altas de forma que f? dz[2m(V — E)]!/2 sea con- 
siderablemente mayor que la unidad) la ecuación (201) es satisfactoria. Sin 
embargo para energías muy pequenas existe otro problema [100]. En el límite 
E — 0 y para potenciales en los que converge la integral iet dz (2mV)?, 
Twp tiende a un límite finito distinto de cero. Por ello el error relativo 
Teract/Twkp puede ser arbitrariamente grande para energías suficientemente 
pequenas. Si queremos ir más allá de la fórmula anterior y obtener la función 
de onda W K B para la colisión con una barrera nos enfrentamos en general 
con el problema de “empate” de soluciones (cinco en un caso típico con dos 
puntos de retorno) debido a la no uniformidad del desarrollo semiclásico (es 
decir a que en cada región se requiere una expresión asintótica diferente) que 
hace a este método un tanto engorroso. Keller [101] propone una representa- 
ción de la solución ui..forme espacialmente que lleva a la siguiente expresión 
para la amplitud de transmisión, 


exp[ió — ia In(a/e)] 
(1 + e272)1/2 


x expi( [^ e — V(z)]? — k)dz 


Io = 
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- (2120) + | {k? - V(z)]"? — k}dz) (202) 


donde 


9. 
II 


argl(1/2+1a), $=0 paraa — 0 (203) 
J "IV (2) — ar. (204) 


ES 
e 
i 


11 Conclusiones 


En esta monografia hemos intentado ofrecer una perspectiva amplia, actua- 
lizada y unificada de las colisiones cuánticas en una dimensión combinando 
resultados conocidos pero dispersos con aportaciones originales. Nuestro de- 
seo es que sea tanto un trabajo de consulta como una incitación a seguir 
profundizando en una serie de problemas abiertos. Confiamos en que el lec- 
tor concluya con nosotros que este no es un campo cerrado y bien establecido 
donde queda poco por hacer. Por el contrario, nos encontramos con multi- 
tud de cuestiones no resueltas de carácter conceptual o técnico relacionadas 
a menudo con aplicaciones de interés tecnológico. La mayoría de estos pro- 
blemas se encuentran en el dominio de estados y/o potenciales dependientes 
del tiempo, que los tratamientos convencionales ignoran o consideran insufi- 
cientemente. 
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